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We suggest a spatially local feedback mechanism for sta- 
bilizing periodic orbits in spatially extended systems. Our 
method, which is based on a comparison between present and 
past states of the system, does not require the external genera- 
tion of an ideal reference state and can suppress both absolute 
and convective instabilities. As an example, we analyze the 
complex Ginzburg-Landau equation in one dimension, show- 
ing how the time-delay feedback enlarges the stability domain 
for travelling waves. 



A common situation encountered in the operation of 
physical systems or devices is that a useful solution of 
the equations of motion turns out to be unstable in a 
parameter regime of interest. In many cases the desired 
behavior is a steady state or a regular periodic motion, 
and the instability eventually leads to chaotic fluctua- 
tions which limit the system's performance. Thus one is 
led to explore the possible modifications of the system 
that render the desired motion stable. 

Recently, there has been intense interest in the appli- 
cation of proportional feedback for stabilizing periodic 
orbits. Since the orbit in question is a solution to the 
equations of motion, the stabilizing feedback signal van- 
ishes when control is successful, so that all the desirable 
features of the uncontrolled system are retained. Many 
methods for controlling systems with only a few relevant 
degrees of freedom have now been successfully demon- 
strated. @|| 

Some of the most interesting and significant dynami- 
cal instabilities arise in spatially extended systems which 
may be described by partial differential equations, a large 
number of coupled ordinary differential equations, or cou- 
pled map lattices. Well-known examples of practical 
interest include convecting fluids, large Fresnel number 
lasers, and arrays of semiconductor lasers. For small sys- 
tems, the number of unstable modes remains small and 
techniques involving only a few degrees of freedom can 
effectively treat the spatiotemporal dynamics. H For the 
case of open systems with convective instabilities, control 
of larger systems has also been demonstrated. Q 

In this paper, we present and analyze a new method for 
stabilizing periodic orbits in arbitrarily large systems. (A 
different approach has been suggested by Hu and Qu. (|]) 
Our approach is a generalization of the technique known 
as "extended time-delay autosynchronization" (ETDAS), 
which has been successfully applied to a variety of low- 



dimensional systems, both numerical and experimental. 
Q In ETDAS, the current state of the system is com- 
pared to its state one or more periods in the past and a 
feedback signal is generated locally (at each spatial point) 
based on the local value of the difference. Our analyt- 
ical treatment of the important special case of the ID 
complex Ginzburg-Landau equation shows both that this 
method can stabilize spatially extended periodic orbits 
and, more generally, that the introduction of spatially lo- 
cal time-delayed interactions can dramatically alter the 
stability properties of extended deterministic systems. 

Our method has three key features: First, it applies 
equally well to systems with absolute or convective in- 
stabilities. Second, it is, by construction, scalable up to 
arbitrarily large system sizes with no increase in complex- 
ity. Third, it does not require comparison to an exter- 
nal reference signal and therefore might be implemented 
in fast (optical) systems, systems in which the reference 
state is not known a priori, or systems in which the ref- 
erence state has nontrivial spatiotemporal structure. In 
systems where the control works, the only information 
that must be supplied by the controller is the period r of 
the desired motion. 

The general approach we take is as follows. To generate 
the feedback signal for a system described by an evolving 
field 4>{x,t), the entire field is compared to time-delayed 
images of itself <j)(x, t — nr), with r chosen to be the pe- 
riod of the desired orbit and n taking all positive integer 
values. With t n = t — nr, the feedback signal is the field 

oo 

E<j,(x, t) = 7 2J R n (<I>(X, t n ) ~ <f>(x, t n+1 )) , (1) 

71=0 

where 7 is a real parameter (the gain) and R is a real 
parameter between and 1. We assume here that r is 
known in advance. (The problem of finding appropriate 
periodic orbits and their periods is beyond the scope of 
this paper.) It is clear that vanishes identically when 
4>{x, t) is periodic in time with period r. As emphasized 
elsewhere Q , the infinite sum can be obtained in practice 
with a recursive feedback loop that contains only a single 
time-delay device. In some systems, it may be possible 
to implement this form of control directly; e.g., by using 
optical elements that preserve the spatial structure of a 
laser beam. Alternatively, e<f,(x,t) can be considered as 
a limiting case of the placement of independent ETDAS 
controllers at many discrete points in the system. 

We are interested in the extent to which proper choices 
of 7 and R can improve the stability of selected time- 
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periodic patterns of the field </>(x, t) for arbitrarily large 
system sizes. In this paper, we treat the ID complex 
Ginzburg-Landau (CGL) equation with a cubic nonlin- 
earity, a partial differential equation that describes a 
large class of systems that undergo a bifurcation from 
regular oscillations to spatiotemporal chaos. [|| In addi- 
tion to its relevance to fluid and laser systems, this equa- 
tion has the advantage of possessing purely sinusoidal 
travelling wave solutions which permit a detailed anal- 
ysis. Our linear stability analysis shows that periodic 
travelling waves states can indeed be stabilized by ap- 
propriate choices of 7 and R, even for system parameters 
corresponding to turbulence ("defect chaos" ||) in the 
uncontrolled equation. An numerical illustration of suc- 
cessful control is shown in Fig. 1. 

The controlled CGL equation we study may be written 
in dimensionless form as 

d t A = A + (1 + ic x )dlA - (1 - ic 3 ) \A\ 2 A + e A , (2) 

where x is a one-dimensional continuous variable, A(x, t) 
is a complex field, c\ and C3 are real parameters, and 
eA(x,t) is the control term defined above. Without the 
control term, Eqn. (Q) admits travelling wave solutions 
of wavenumber k and frequency lu = — C3 + (ci + c 3 )k 2 . 
Each solution 



the Floquet multipliers of a general, finite-dimensional 
system controlled by ETDAS is, 01 



Ak(x,t) = yl — k 2 exp(ikx — iut) 



(3) 



becomes unstable for large enough c\ and/or C3, and all 
of them are unstable for C1C3 > 1. ||] 

We find that when r is chosen to be 2-k/lu, the domain 
of ci and C3 values over which the solution Ak is stable 
is expanded significantly for modest values of 7 and R. 
Some typical results are shown in Fig. 2. Each panel 
of the figure corresponds to a different set of values for 
k and R. The shaded region represents the parameter 
values for which the travelling wave Ak is a stable solu- 
tion of the equation with control, but would be unstable 
without control. The instability would be convective just 
above onset, but absolute for larger c\ or C3 [Q, as indi- 
cated by the dashed line in Fig. 2. A surprising result is 
that time-delay feedback allows stabilization of travelling 
waves deep into the ordinarily chaotic, absolutely unsta- 
ble regime [c\c 3 ~ 8), even though it is almost totally 
ineffective in stabilizing the uniform (k = 0) oscillatory 
state. 

We now describe our procedure for obtaining the sta- 
bility domains depicted in Fig. 2. Standard linearization 
of Eq. (||) about Ak{x,t) yields sets of ordinary, time- 
delay differential equations for the Fourier amplitudes of 
a perturbation. The technique of Ref. Pj is then applied 
to determine the stability of the different modes. In each 
period t, a given mode grows or decays by a complex 
factor n (a Floquet multiplier) . A system is stable if and 
only if < 1 for every mode. The defining relation for 



X dt ( J+7 - 
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0. 
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where T[. . .] represents the time-ordered product, J is 
the Jacobian of the uncontrolled mode equations, M is a 
"control matrix" that contains the information about the 
way in which the control signal is formed and enters into 
the dynamical equations, and 1 is the identity matrix. 

In the present case, three features simplify the analy- 
sis: first, M = 1; second, the Fourier modes decouple, 
with each yielding a condition of the form of Eqn. (jij) 
with 2x2 matrices J and M; third, neither J nor M is 
time dependent. The latter is due first to the trivial time- 
dependence of the desired solution Ak and second to the 
directly additive way in which e appears in the equation 
for the controlled system. In such cases, the determinant 
in Eqn. (^) may be evaluated explicitly by solving the 
differential equation that yields the time-ordered prod- 
uct. Here the defining relation for the Floquet multipli- 
ers associated with a perturbation at wavenumber q + k 
becomes, 

gifi' 1 ) = fi- 2 e 2aT - 2 t r 1 e aT cosh (/3r) + 1 = 0; (5) 



2icikq + k 2 - 1 + 7 



1-Rfi- 1 ' 

13 = [(1 - k 2 - Aic 3 kq + 2 Cl c 3 q 2 ) (l - k 2 ) 
+4k 2 q 2 + 4i Cl kq 3 - C?Q 4 ]5. 

The state Ak(x,t) with particular choices of c±, C3, 7, 
and R is linearly stable if and only if all of the roots of g 
lie outside the unit circle. Note that has an infi- 

nite number of roots due to the time delay in the system. 
As described in Ref. || > it is straightforward to perform 
a winding number calculation (or a contour integration) 
that will return A/", the number of roots that lie inside the 
unit circle. The linear stability condition then reduces to 
M = 0. The winding number calculation is performed nu- 
merically by evaluating g on selected points on the unit 
circle. The precision of such a technique is determined by 
how well one can distinguish between a root inside the 
unit circle and one very close, but still outside. Using 
an adaptive step-size method, we resolved the location of 
such roots to an accuracy of 10~ 6 . 

It is at this point that the many degrees of freedom 
in a spatially extended system complicate the analysis. 
In order for a particular state to be stable, it must be 
true that a single value of 7 exists for which Ak is stable 
with respect to perturbations at all wavenumbers. To 
see whether such a 7 exists for fixed k, c±, C3, and R, it 
is useful to plot the region of stability in the space of 7 
and the perturbation wavenumber q. Jl0[ Fig. 3(a) shows 
an example for which Ak is linearly stable against per- 
turbations of all wavenumbers for a range of 7 (shown 
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between the dashed lines). Note that the plot must be 
symmetric about q — since from Eqn. || it is clear that 
n(q) — fx*(—q). The rapid divergence of the stability 
boundaries for large q merely reflects the fact that the 
system is highly stable with respect to large q perturba- 
tions in the absence of control. 

Figs. 3(b) and (c) show why control cannot be achieved 
for some values of c\, C3, k, and R. The problem is that 
peaks in the lower boundary reach to values of 7 that are 
already ruled out by valleys in the upper boundary, so 
that no single value of 7 can stabilize all wavenumbers. 
The source of the peaks may be understood as follows: 
For a periodic state with frequency w and J and M in- 
dependent of time, it can be shown that ETDAS cannot 
stabilize a perturbation for which arg /1 = mw, where m 
is any integer. The peaks in the lower boundary occur at 
wavenumbers where this condition is approximately sat- 
isfied. In Fig. 3, the peak in the lower boundary at q = 
corresponds to m = 0. In Fig. 3(c), the peak at q ~ 0.75 
corresponds to m = ±1. In the present case there can 
be no effect from higher |m| because control will already 
have been lost due to m = ±1. 

By analyzing stability diagrams in the 9 — 7 plane for 
a grid of values in the c\ — C3 plane, one can construct 
the stability diagrams shown in Fig. 2. As in the ex- 
ample of Fig. 3(a), there is a range of the feedback gain 
7 that successfully achieve control for each point in the 
shaded area of Fig. 2. In general, the minimum value of 
I7I required increases smoothly as c\ and/or C3 increases. 
In the domains shown here, as well as others we investi- 
gated, the value of I7I required for stability is less than 
1 even at the highest values of c\ and C3 in the control- 
lable domain. (Details will be given in a longer paper.) 
Beyond the line labelled U0 (Ul) control is lost due to 
the mechanism described above with m = (m — ±1), 
not through a divergence in the required 7. 

Fig. 2(a), (c), and (d) illustrate how the stability 
boundaries shift for different choices of k. Fig. 2(a) cor- 
responds to k = .075-7T. For larger k (Fig. 2(c)), the U0 
line moves farther from the uncontrolled stability bound- 
ary and Ul moves closer. For smaller k (Fig. 2(d)), the 
situation is reversed. As k is decreased toward 0, the 
boundary U0 approaches the original uncontrolled stabil- 
ity boundary, so that no enhacement of the k = state is 
obtainable. Figs. 2(a) and (b) show the effect of changing 
R. As R is increased, the domain of stability increases in 
area. However, even as R approaches its maximum value 
of 1, the domain of stability cannot include the region 
in which one of the unstable modes of the uncontrolled 
system has frequency ±w. The boundary of this region 
is the dotted line in Fig. 2(a). 

We have checked specific aspects of the results pre- 
sented in Fig. 2 with numerical simulations of the con- 
trolled CGL equation. Periodic boundary conditions 
were employed with the system size chosen to be an in- 
teger multiple of the wavelength of the travelling wave. 



System sizes corresponded to a length of at least 15 x 
2ir/k. The simulations were performed with a second- 
order predictor-corrector and finite difference technique, 
with time steps of order 10~ 2 and spatial resolution ~ 400 
points. The instabilities were observed to occur at values 
of (01,03) consistent with the analytic results presented 
here for infinite systems. 

We have demonstrated that time-delay feedback can 
be effective in stabilizing periodic states of spatially ex- 
tended systems. Application of this technique to the sta- 
bilization of unstable ordered states in fluid, laser, and 
biological systems is strongly suggested. We expect the 
control technique to be applicable to many types of pe- 
riodic states, though the stability analysis may become 
complicated. If the linearized equation for the pertur- 
bations about the periodic solution has space-dependent 
coefficients but its time dependence is still trivial, per- 
turbations can be decomposed into appropriate eigen- 
functions and the analysis discussed here will apply. If 
the periodic state has trivial spatial dependence but non- 
trivial time dependence, then the stability of the Fourier 
modes can be analyzed using the numerical method of 
Ref. H . Finally, when the periodic state has complicated 
spatiotemporal structure, it appears that numerical inte- 
gration of the controlled equations would be the most 
efficient approach. Even in the absence of any stability 
analysis, however, control can be attempted in a physical 
system given only the knowledge of r and the ability to 
adjust the single parameter 7. 

Our work points to several important questions for fu- 
ture study. What is the minimum density of discrete 
controllers needed in situations where spatially continu- 
ous processing in the feedback loop is not possible? What 
level of noise can be tolerated? How can one force the 
system from the spatiotemporally chaotic state into the 
desired controllable state? 

Finally, we suggest that the application of time-delayed 
feedback may be a valuable tool for studying the intrinsic 
physics of spatiotemporally chaotic systems. By varying 
7 and t slowly, it may be possible to locate previously 
unknown periodic states or to observe other novel effects. 
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Figure 1. Numerical simulation of the ID CGL equation 
with periodic boundaries and parameters (ci = C3 = 2). 
(See final page of text for simulation details.) The plots 
show the evolution of the phase of the complex field A. 
Black (white) regions correspond to phases near (±7r). 

(a) Without control, an unstable travelling wave state 
(wavenumber k = .0757r) evolves into a turbulent state. 

(b) Time-delayed control stabilizes the travelling wave. 
The control parameters (see Eq. (1)) are 7 = —.6 and 
R = .75. 



Figure 2. Stability diagrams for three choices of wavenum- 
ber k and two choices of parameter R in Eqn. (1). (a) 
k = .0757T, R = .75, The circles mark the values of c\, C3 
used in Fig. 3. (See final page of text for explanations 
of UO, Ul, and the dotted line.) (b) k = .075tt, R = .5, 
(c) k = .1tt,R = .75, (d) k = .05ir,R = .75. In the area 
labelled S the uncontrolled travelling wave is stable. In 
the shaded region, control is possible. In all cases shown 
here, control is achieved with gain I7I < 1. In the re- 
gion marked U, no value of 7 can stabilize the travelling 
wave. The dashed line marks the transition from con- 
vective (lower left region) to absolute instability in the 
uncontrolled system. 
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Figure 3. Maps of the stable domains (shaded area) for 
k = .0757T and R = .75 in terms of the feedback gain 7 and 
the perturbation wavenumber q for four values of (ci, C3). 
(a) ci = C3 = 2 (b) ci = C3 = 2.5 (c) ci = 3.3, C3 = 1.75. 
The dotted lines are drawn at 7 = to draw attention to 
the case of no control. In (a) the dashed lines show the 
range of 7 for which control is successful. 



1 



